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Abstract 



We give a precise numerical solution for decorated Ising models on the simple cubic 
lattice which show ferromagnetism, compensation points, and reentrant behaviour, 
ly-j . The models, consisting of S = \ spins on a simple cubic lattice, and decorating 

Q\ [ S = 1 or S = | spins on the bonds, can be mapped exactly onto the normal spin-^ 

' Ising model, whose properties are well known. 
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Ferrimagnets are materials where ions on different sublattices have opposing 
magnetic moments which do not exactly cancel even at zero temperature. [1] 
An intriguing possibility then is the existence of a compensation point, below 
the Curie temperature, where the net moment changes sign. This has obvious 
technological significance. 

There has been considerable work in recent years in studying these phenom- 
ena through simple models[2,3,4,5,6,7], where treatments beyond mean-field 
theory are possible. Of particular interest are decorated systems, which can 
be mapped exactly onto simpler models, and in this way solved either exactly 
or to a high degree of numerical precision. 

The study of decorated Ising model or, more generally, of Ising model trans- 
formations has a long history [8, 9]. Kaneyoshi[5] introduced a mixed spin Ising 
model on the square lattice with spins Sa = \ on the vertices, coupled antifer- 
romagnetically to spins Sb > \ decorating the bonds. Within the framework 
of an effective-field theory with correlations, a number of interesting results 
were obtained. In particular, it was found that for Sb = 1 and a negative crys- 
tal field anisotropy term the system could show up to three separate phase 
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transitions and one or two compensation points. The crystal field term and 
a residual ferromagnetic coupling between the A spins were crucial for these 
effects. On the other hand there was no reentrant behaviour and at most one 
compensation point for Sb — §■ Dakhama[6] obtained an exact solution for 
this same system, and confirmed the existence of compensation points, as pre- 
dicted by Kaneyoshi[5], but found no evidence for multiple phase transitions. 




• Spin \ — J 
O Spin S --- J' 

Fig. 1. The mixed spin cubic lattice. 

It is quite straightforward to extend this analysis to the simple cubic lattice, 
and this is the substance of our paper. Indeed it seems surprising that this 
has not been done before, but to our knowledge it has not. We consider a 
simple cubic lattice with N spins Sa = \ on the vertices and 3N spins S B > \ 
decorating the bonds. The model and its interactions are shown in Figure 1. 
The Hamiltonian is taken as 

H = J ]T cnSj -J'Y, Wi - A £ S? (1) 

(ij) (ki) i 

where a { = ±\, S { = -S, -S + 1, • • • , S and J, J' > 0. 

Summing over the B spins we can write the Boltzmann weight for each bond 
as 

s 

exp[if'(7i(72 - KfJlfa + (7 2 ) + Dfl 2 ] = A exp(K<7i<7 2 ) (2) 

fi=-S 



where we have introduced coupling constants K = f3J, K' = /3J', D = j3A. 
Thus the model transforms to a spin-| Ising ferro magnet on a simple cubic 
lattice, with coupling constant 

K = ^' + 21n(S 1 /S 2 ) (3) 



2 



with 



D/i 



(4) 



The prefactor A, which is needed to compute the free energy, is given by 
A 2 = £iE 2 . Note that the new coupling constant K is independent of the 
sign of J, and hence our conclusions regarding phase transitions hold equally 
well for ferromagnetic nearest neighbour coupling. Of course, in that case, the 
sublattice moments will be parallel and there will be no compensation point. 

The critical point of the nearest neighbour spin-| Ising model on the simple 
cubic lattice is known to high accuracy [10] as 

K c = 0.88664 (5) 



(note that we use a — ±|, rather than a — ±1) and hence our model has a 
critical line, which can be expressed as 



e K'/2 



]T e D " 2 cosh Kfi = 1.5579 £ e Di? (6) 



For the two cases which we consider explicitly, namely S — 1, § we can write 
this is 

S = l: e K ' /2 {l + 2e D cosh K) = 1.5579(1 + 2e D ) (7) 

S=-: e^cosh - + e 2D cosh — ) = 1.5579(1 + e 2D ) (8) 
2 2 2 

These equations must be solved numerically, and the solutions are discussed 
below. We note, of course, that Eqns. (6)- (8) have the same form for any 
lattice, with the appropriate e Kc ^ 2 . 

The spontaneous magnetizations of the two sublattices are easily obtained. 
For the A sublattice we have simply 

(a) = l -M {K) (9) 



where Mq(K) is the magnetization of the spin-| Ising model on the simple 
cubic lattice, which can be obtained numerically to high accuracy. The mag- 
netization of the B sublattice is given by (precisely as in [6]) 

(S) = 2 7 < a > (10) 
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where 

7 = £ 3 /£i (11) 



and 



£ 3 = J2ve D ^ sinhif/i 



(12) 



The total magnetization is then given by 



M/N=(a) -3(S) 

= (l-67)<(7> (13) 

whence any compensation point must lie on the line 67 = 1. Note that this 
line depends on D but not on the second neighbour interaction K'. Explicitly 
for S — 1, § the compensation point line is 



S=l: i2e D sinhfr = 1 + 2e D coshK 



3 n ■ 1 % 

- : 6 sinn — + 18e 



2D 



sinh 



3K 



cosh — + e 



2D 



cosh 



(14) 
(15) 



We now turn to the results, which we stress are exact to within the very small 
numerical uncertainty in K c . 

We first consider the case 5 = 1. Figure 2 shows critical lines in the K^ 1 — 5 
plane (5 = D/K) for various values of a — K'/K. We see, as expected, that 
increasing the second-neighbor coupling increases the critical temperature as 
does increasing the anisotropy parameter A. The critical line for a = shows 
reentrant behaviour, i.e. two phase transitions, for a small range —1.14046 < 
5 < —1.0. At very low temperatures the system will be disordered, since 
A < favors that state S = at B sites. As the temperature increases the 
system orders at a lower critical temperature T^, passes through an ordered 
phase, and then disorders again at an upper critical temperature Tc 2 . Any 
nonzero second neighbor interaction will give rise ordering at T = and at 
low temperatures, and hence there will be a single critical point (subject to the 
caveat below). Reentrant behaviour of this kind was found in the approximate 
treatment of the square lattice[5], but was shown to be spurious by the exact 
solution [6]. For the simple cubic lattice it is a real effect. 

The compensation point line is shown in Figure 2 as a dotted line. If this line 
lies below the phase transition line a compensation point will be present. Thus 
for a = these is no compensation point, i.e. a second-neighbour interaction 
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Fig. 2. The phase diagram of the mixed spin cubic lattice for Sb = 1- Lines give 
critical temperature versus single ion anisotropy for various strengths a of the sec- 
ond-neighbour interaction. The dotted line is the compensation point line (Eq. 14). 

is essential. For 5 < —2.7858 no compensation point is possible, while for 
—2.7858 < 5 < — 1 two compensation points will occur for sufficiently large a. 
For 5 > — 1 a single compensation point will occur for sufficiently large a. 
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Fig. 3. Critical temperature versus anisotropy for small a. In the range 
< a < 0.1152 the system will have three critical points. The dotted line is the 
compensation point line (Eq. 14). 

Careful analysis of eqn.(7) for small at reveals even more interesting behaviour. 
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In Figure 3 we show critical temperature curves in the vicinity of 5 = — 1 for 
several a values. As is seen, for < a < 0.1152 the curve is S-shaped, and 
hence the system has three critical points, going from an ordered state to 
a disordered state, then back to an ordered state and eventually to a high- 
temperature disordered state. While such behaviour has been seen before in 
approximate treatments, as for example in the square lattice by Kaneyoshi[5], 
we are unware of any cases where it has been demonstrated in an exact treat- 
ment. 




Fig. 4. (a)Total magnetization versus temperature for a = 3.2, 5 = —2.5, showing 
two compensation points; (b)Total magnetization versus temperature for a = 0.05, 
5 = —1.2, showing three critical points. 

It is useful to compute the total magnetization as a function of temperature, 
for particular (a, 5), to display some of the features discussed above. The 
magnetization M (K) for the simple cubic lattice is not known exactly, but can 
be computed numerically to high accuracy using the existing low-temperature 
series[ll] and the best estimate of the critical behaviour [12]. In Figure 4(a) we 
shows magnetization versus temperature for the case a = 3.2, 5 = —2.5. The 
two compensation points are evident. Figure 4(b) shows magnetization versus 
temperature for a = 0.05, 5 = —1.2. This case has three critical points. The 
compensation point line falls in the intermediate disordered phase, and hence 
there is no compensation point for this case. There are parameter values where 
the compensation point line falls in the lower T ordered phase, giving three 
critical points and a compensation point: an example of this is for a = 0.11, 
5 = —1.3343, where we have three critical points at K = 3.12055, 4.18638 and 
4.60271, respectively, while the compensation point is at K = 4.81408, inside 
the lower T ordered phase. 

Figure 5 shows the critical point curves and compensation point curve for 
s — \. There is a clear qualitative difference between the two cases S — 1 and 
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S — §, as previously seen for the square lattice [6]. In the S — | case there is no 
compensation point for a < 2.5519, for any 8, whereas a single compensation 
point occurs for any a > 2.5519, for a range of 5 values. In all cases there is a 
single phase transition point. 

In conclusion, we have extended the analysis of Ref . [6] to the case of mixed spin 
Ising models on a decorated simple cubic lattice. Using the numerical precise 
value of the critical point for this lattice we obtain essentially exact results for 
phase transitions and compensation points. For the S — (|, 1) system we find 
reentrant behaviour, with the possibility of three, two and one phase transition 
and zero, one, or two compensation points. For the S = (|, |) system the 
behaviour is qualitatively different, with always a single transition and zero 
or one compensation point. This difference can be qualitatively understood 
through population of the S = state at low temperatures, and we expect it 
to persist for higher integer and half-integer spin values. 

For the models studied here, the presence of second-neighbour interaction is 
essential for the occurrence of a compensation point. The same conclusion, for 
a different mixed spin Ising model, was noted in Ref. [4]. 

It is interesting to note that our (essentially) exact results are qualitatively 
similar to those obtained by a correlated effective field theory for the square 
lattice[5]. Presumably the reduced fluctuations in 3-dimensions make the effective- 
field theory more reliable than it is for the square lattice. 
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Fig. 5. The phase diagram for Sb = §■ Lines give critical temperature versus single 
ion anisotropy for various strengths a of the second-neighbour interaction. The 
dotted line is the compensation point line (Eq. 15). 
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